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We construct a theory of quantum fluctuatons in a regular array of small superconductive islands 
of size d connected via low-resistance tunnel contacts (Gt = h/4e 2 R t 3> 1) to dirty thin metal film 
with dimensionless conductance j > 1. Electrons in the film interact repulsively with dimensionless 
strength A. The system is macroscopically superconductive when distance b between neighbouring 
islands is short enough. The zero-temperature phase transition from superconductive to normal- 
conductive state is shown to occur with the increase of distance between superconductive islands, 
at Inbc/d ~ G\j\g. The critical distance b c is much less than 2d localization length L; oc ~ e 19 , 
so the considered effect develops when weak- localization corrections are still small. The T c (g,b) 
dependence at b < b c is found. These results are valid at sufficiently large g, whereas decrease of g 
is expected to lead eventually to another b c (g) dependence, lnfe c /d ~ yfg. 



1. Introduction. Quantum phase transitions take place at zero temperature due to variation of the parameters of 
the system. An important group of examples of such phase transitions is related to superconductors. The character of 
a phase transition here depends on the type of the system. In granular superconductors, a phase transition between 
macroscopically superconductive and insulating states occurs when the Josephson coupling energy becomes less than 
the Coulomb energy due to transfer of a single Cooper pair between two grains [Q . Experiments demonstrating this 
type of phase transition in thin superconductive films are numerous . More recently, similar results were obtained 
on artificially prepared almost regular Josephson arrays (^-j^) . A theory of superconductor-insulator phase transition 
was discussed in HQ. A different kind of 2D superconductor - homogeneously disordered films - was considered by 
Finkelstein Q who predicted superconductive-normal phase transition (vanishing of T c ) due to disorder-enchanced 
Coulomb repulsion in the Cooper channel (relevant experimental results are discussed in ||). 

Within BCS approximation nonmagnetic impurities do not suppress the temperature of superconductive transition 
("Anderson theorem", cf. []To| , [ll|| ) . However in the situation when T c is suppressed by quantum fluctuations, non- 
magnetic impurities can lead to further suppression of T c . An example of such an effect in quasi-one-dimensional 
superconductors was considered by V.I.Mel'nikov and one of the authors Another example of this kind is studied 
in the present paper. We consider 2D superconductive system where small superconductive (SC) islands are con- 
nected via insulating barriers to a dirty metal film with good conductance, g — a ■ 27rS,/4e 2 3> 1 (below we put fi = 1). 
Normal-state tunnel conductance a t — Gt • 4e 2 /27r?i of the barrier is assumed to be high, Gt 3> 1. We show that 
in this case quantum fluctuations destroy superconductive order even at T = 0, if distance between neighbouring 
SC islands b becomes sufficiently long (in this paper we will assume regular lattice of SC islands). Although the 
critical inter- island distance b c appears to be exponentially long at g ^> 1, Gt 3> 1, it is still much shorter than the 
2D localization length L/ oc ~ e~xjp(ng), so the destruction of superconductivity produces normal-metal state (at least 
above weak localization temperature T/ oc ~ D/L 2 oc ). 

Since we consider temperatures much below local T® of SC islands, any fluctuations of their gap magnitude |A| can 
be neglected, and the "phase-only" Hamiltonian is appropriate, with its part responsible for the inter-island coupling 
due to proximity effect in the normal conductor (or due to direct tunnelling as in Josephson junction arrays) given by 

Hint = ~ ^ Jjj cos(0. t - <pj) (1) 

ij 

If the coupling goes via 2D metal, then at T = and neglecting electron-electron interaction, one gets (cf. e.g. [[r3[) 
Jij = J{\rij\) tx l/\rij\ 2 . Such a very slow decay of coupling strength allows us to treat the inter-island coupling in the 
mean- field approximation (MFA). If the temperature is not sufficiently low, then only nearest-neighbouring islands 
interact and MFA can be used for an order-of-magntiude estimates of the transition point. Within MFA the action 
of all other islands on a given one is described by the mean "field" they produce: J(cos^) J dt cos </)(t). As a result, 
the transition point for the appearance of superconductivity at low temperatures is determined by the condition (cf. 
e.g. where a similar transition in quasi-lD system was considered) 
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1 = \j f ' C (t)dt; J = J2 Jij ; C (t) = (cos(0(O) - <f>(t))) = 2(cos0(O)cos0(i)) (2) 
1 J ° i 

where Co(t) is the imaginary-time autocorrelation function describing phase fluctuations of a single SC island in the 
metal film. These phase fluctuations are determined by the charge transport between the SC island and surrounding 2D 
metal. Since we are going to consider very low temperatures T« A and, moreover, the energy scale of all phenomena 
related to establishing the inter-island SC coherence is much below A, the only channel of charge transport between 
SC island and 2D metal is the Andreev subgap conduction | jl5| . Corresponding phase fluctuations can be conviniently 
described by the imaginary-time dissipative action introduced in jl6|,[l^] , which is similar to the one of single-electron 
tunnelling in normal tunnelling junctions (cf. e.g. |l8|]): 

SoW)] =-J dtdt'K A {t - t') cosM*) - 0(f)); K A (t) = J^-j (3) 

where Ga is the Andreev subgap conductance (normalized to conductance quantum 4e 2 /27r?i) between SC island and 
2D metal. We will always assume that Ga S> 1. 

The crucial feature of this dissipative action as compared to the standard dissipative action for Josephson junctions 
with resistive shunt |l9[] , is that it is periodic in the phase difference <fi(t) — <fi(t'), thus it correctly describes the fact 
that the electron charge in the island can be changed by 2e quanta only. As a result, in the longest-time limit t t c 
the correlation function Co(t) decays sufficiently fast (~ 1/i 2 ) and the integral J °° Co{t)dt oc t c always converges, 
which is crucial (cf. Eq.(||) for the existence of the zero-T quantum phase transition. The "correlation time" t c 
introduced above appears to be exponentially long, t c oc e° A l 2 (cf. §^,|l|). It means that the Josephson coupling 
must be exponentially weak in Ga in order to destroy superconductive coherence in the array at T = 0. On the 
contrary, in the linear dissipation model JToj ] coherence is preserved at arbitrarily weak J as long as the dissipative 
conductance Gat > 1 (at large Gn the correlation function decays very slowly, C (t) ~ i- 2 /Giv ). 

Below we will analize (in Sec. 2) the behaviour of Andreev conductance Ga as a function of the normal-state barrier 
conductance Gt and the sheet conductance of the metal g. Then we will calculate (in Sec. 3) the low-temperature 
phase boundary between macroscopically superconductive and normal-metal states of the array. The Sec. 4 is devoted 
to the discussion of our results and to the conclusions. 

2. Andreev conductance and phase correlations of a single island. 

We start from the case of relatively low transparance of the insulating barrier between SC island and 2D metal film, 
i.e., we will assume that the subgap Andreev conductance Ga can be calculated to the second (i.e. lowest nonzero) 
order in the normal-state tunnelling probability (presise condition will be discussed below). Convenient expression 
for the phase-coherent subgap conductance can be found, e.g., in Eq. (12) of paper [^2| : 

GUM = ^ln^ (4) 

where Gt = 2-kTi j (2e) 2 Rt is the dimensionless tunnelling conductance in the normal state, uid = D/d 2 , with d being 
size of the SC island, and uj is the frequency of the measurement (for dc conductance, uj should be replaced in Eq.(|J) 
by the applied voltage V). Extra (with respect to ||^]) coefficient 2 in the denominator of Eq.(||) is due to a different 
geometry: in our case Cooperon propagates over the whole 2D metallic plane, whereas it was half-plane in Eq.(12) 
of p2|| . The physical meaning of the factor {l/g)hi^ in the expression (jij) for Ga is that it provides the return 
probability for two electrons to come back to the same point where they tunnel from the SC island, i.e. it is the time 
integral over the Cooperon return amplitude C(r = 0,t). Below we will always assume that Ga{u ~ <-^d) 3* 1, i-c 

G° A = ^»l (5) 
2ng 

Expression (||) was derived assuming no electron-electron interaction in the normal metal. It can be generalized to 
the case of nonzero (dimensionless) repulsion in the Cooper channel A > 0. To do so, we note following [|l3| that in 
the presence of interaction 

1 ln^ 

dte lwt C(r = 0,t) = —- (6) 

2ng 1 + A In w 

where 
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A = A(d) = ^— ^- (7) 

1 + A ln^ 

is an effective dimensionless constant of electron-electron repulsion in the metal (denned on the scale d) . Therefore 
at sufficiently low frequencies u> < oj\ — ujde^ 1 ^ the subgap conductance becomes frequency-independent, 

G A = ^- y *<o* (8) 

The results (||j^) for subgap conductivity were derived under the assumption that the superconductive phase </> of 
the island does not fluctuate ( due to the coupling to external low-impedance circuits). In our problem we deal with 
an array of many small SC islands which are coupled to the surrounding only by the Andreev conduction processes, 
therefore we should take phase fluctuations into account explicitely. We will calculate, using the dissipative action of 
the type of (||), the single- island correlation function Co{t) as defined in Eq.(||) and the effective charging energy Eg 
defined via the relation 

/>oo 

1/E* C = / C (t)dt (9) 
Jo 

We start from the assumption that the interaction constant A is not very small, so the relevant frequency scale 
u> <C and Andreev conductance Ga is frequency- independent and given by Eq.(||) The corresponding action 
is given by Eq.(3h, which is formally equivalent to the problem considered by Kosterlitz p3| and then later on by 
others |^0|,[l8 24 2[| . The bare correlation function for the phase fluctuations as determined by Eq. (|J) is given by 



w»' 2 > = wr A < 10 > 

and leads to the logarithmic growth of phase fluctuations in the time domain at t 3> cl^ 1 (upper limit of this region 
will be discussed below), 

So(t) = \((Ht) - 0(O)) 2 )o = (2/GU) \nu x t (11) 

Now one could calculate Co(t) in Gaussian approximation using Eq.(|lT|) and get power-law decay with the exponent 
2/Ga', however, such a procedure would be incorrect since the parameter Ga gets renormalizcd due to the nonlinearity 
of the action (|^). The renormalization group (RG) equations to the main order in 1/Ga were derived in |p3| : 

dG A (Q _ d\nCo(C) _ 2_ 

d( ' d( G A (0 [ ' 

where £ = m). These equations can be trivially integrated, the solution is 

Ga (w) = Ga — 2 In — ; C (t) = 1 - -|- lnw A i (13) 
uj Ga 

Physical meaning of Ga(w) dependence as given by Eq.(|l^) can be understood as follows: due to the fluctuations 
of 4>{t), the coherence between successive Andreev reflections of the electron which have been diffusing in the metal 
during long time ~ 1/w, is partially destroyed, which leads to the decrease of subgap conductance on the corresponding 
timescale. Note that the Eqs.( p^Jl^) are valid as long as renormalized Andreev conductance is large, Ga{C) ^ 1) i- e - 
Co(t) decreases according to Eq.fll3j) till t » t c = co^ 1 e GA ^ 2 , where Co(t) becomes of order 1/Ga <C 1. At longer 
timescales t 3> t c the correlation function Co(t) ~ G^ 1 (t c /t) 2 . The main contribution to the integral J °° Co(t)dt = 
l/Eg comes from t ~ t c . Using ( |l3| ) we obtain the estimate E* c « ^ ss u>de~ 1 ^ x GAe~ GA ^ 2 . The above evaluation of 
t c was based on the leading order of RG approach, thus it does not, in general, give the correct preexponential factor. 
A more accurate calculation of E% in a similar problem was done in with the use of the instanton technique 
developed in |21|. The result of (confirmed later on in the Monte Carlo simulation |2^]) can be translated in our 



case as 



E* c = A» d e-^ 2 e- G -/ 2 = Alo^ (0- ) e -&V^ (14) 
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where A is some numerical factor of order 1. Thus the effective Coulomb energy E* c associated with the Andreev 
conduction process is still nonzero, although exponentially suppressed at large G a (basically the same result was 
mentioned in |l7j]). This is why we expect quantum ( at T = 0) phase transition to occur when the Josephson 
coupling energy J is of the order of Eq. 

Before we turn to the discussion of the phase transition, let us analize phase fluctuations in the frequency regions 
where bare Andreev conductance differs from (JsJ) , and is frequency-dependent (cf. Eqs.(|4|). The corresponding 
dissipative action is of the form (||), but with a generalized expression for the kernel K A (t), which is defined via 

J dtK A (t)(e 1 ^ - 1) = -Mg a («) (15) 

Now, Eq.@ will be the same but with frequency-dependent Ga(u>), which leads to the following results for So(t): 

S o(t) = -^Inlnwdi; oj^ 1 < t < oj^ 1 (16) 

In the shortest time window t < w, 1 fluctuations grow extremely slowly, and can be neglected altogether due to the 
smallness of the prefactor (cf. Eq.(|5])). 

3. Coupling due to proximity effect and the SC-M transition. 

Total action of the system contains both single-island dissipative terms determining phase dynamics of individual 
islands, and Josephson coupling between phases on different islands 



S = - I dtdt'Y^ - *0 cos(0j(t) - ^(t 7 )) - / dtdt'^2 J jk (t- Ocos(^(t) - 4> k {t')) 



(17) 



As was explained in the Introduction, we are going to look for the SC-M phase transition within the mean field 
approximation, which is possible due to slow decay of the Josephson coupling energy J Jjk(t)dt = J{Tjk) with island 
separation r^- Josephson energy can be found (in the leading order in G t ) with the use of the results |li| (note that 
in that paper planar geometry of S-N-S contact was considered, whereas now we deal with the cylindrical geometry). 
Then we come to the following expression for the Fourrier component of the coupling 3(<j): 

J{q) = ^YTWqV ^ = 2 " T ^ 2Ki + W = ln max^D?) (18) 

Note that neglecting electron-electron interaction A one would get logarithmic divergence of J{Q) in the zero-T limit. 
The expression (|l^) for the Josephson coupling is analogous to Eq.(||) for the subgap conductivity in the presence of 
interaction A. We will always assume Aln ^ 3> 1. In the real space, Josephson interaction is given (at T = 0) by 

G 2 1 

JW = 27rW(l + 2Aln^) 2 (19) 

We see that J(r) does not depend on the metal conductance g. We need to find the sum J = J( r n) ov er the 
lattice of islands. This sum can be approximated by the integral over d 2 r with the lower limit r m i n ~ b: 

J = L — J{r) = ^¥T^ d (20) 

Using now Eqs.(j|,^,|lJ,^), we obtain the critical value of inter- island distance b c such that at b > b c macroscopically 
superconductive state is absent even at T = 0: 



b c 7r 2itg^/~\ 
~d = V2A G 



2 — exp 
t 



2A \ing J 



(21) 



The above formula for b c may look a bit paradoxical, since it shows that the critical distance between SC islands 
grows with the decrease of the film conductance g. Of course, such a result is valid in the intermediate range of not 
very small g, where both Eq and J can be calculated in the leading order in Gt- Now we discuss the conditions for 
such a procedure to be valid. 
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An important condition for the validity of Eqs.([l3 14) comes from consideration of the low- frequency region. First 
of all, the expressions (|^J|) are written in the lowest order in G t and thus the necessary condition for their validity 
consists in the smallness of the anomalous quasiclassical Green function, F(r, oj) <C 1 everythere in the normal metal, 
including the vicinity of the tunnel barrier r = d. If this condition (which is equivalent to the inequality G a *C Gt) 
is not fulfilled, the simple perturbative theory in terms of tunneling amplitude becomes invalid, and a full nonlinear 
treatment based on the Usadel equations for the Green functions F(r, to) and G(r, to) must be developed (see e.g. Q 
and references therein). Thus the necessary condition can be written as 

which must be fulfilled together with Eq.(||). 

There is also another assumption implicitely used in the derivation of Eq.(|2l]). Namely, we neglected Coulomb- 
induced contribution || to the Cooperon-channel coupling constant A, which is of the order of S\ ~ (l/<?)ln^. 
The estimate for 5X may be derived also along the lines of Ref. 27, where Coulomb correction 5H to the Cooperon 
amplitude was calculated in the main order over 1/g. The results of j27j, considered in the limit T — > 0, shows that 
the relative correction 

— — In — . (23) 

11(cj) g uj 

Comparison between Eq.(^3|) and Eq.(^J) shows again that Coulomb corrections can be neglected as long as 

A»iln^ (24) 
g to 

Relevant scale of frequency is given by i" 1 = u\e~ GA ^ 2 , so the estimate Eq.(p4[) actually coinsides with the condition 
given by Eq. (p2|) . Therefore the inequility ( p2| ) is needed for both our assumptions (weakness of tunnelling and of 
Coulomb interaction in the metal) to be valid. 

Currently we can only speculate about the opposite case A < Gt/2irg, which is necessarily realized as g decreases. 
It is clear from the previous estimates, that in this case both Coulomb interaction in the metal and higher orders in 
tunnelling (usually handled by means of the Usadel equation) should be taken into acocunt. For the purely quialitative 
discussion, let us first put A = and consider the results [^6| for the subgap conductance beyond the perturbation 
expansion in Gt- According to po| , the behaviour of the total subgap conductance Ga(<-o) is determined by the 
product G t r 2 D, where t 2 d is the total resistance of the diffusive part of the conducting circuit. In our case it is the 
resistance of the metal film measured between the SC island and the second contact assumed to be big circle of radius 
R > d. Thus r 2D = (l/2ng) In L{u J R) , there L(lu,R) = min(^£>/w, R)- If G t r 2D < 1, we are back to the result 
given by Eq.(Q). In the opposite limit Gtr 2 o ^S> 1 the total subgap conductance was found to be determined by the 
conductance of the metal film only and is insensitive to the barrier transparency: 

GaM = r4£b), (25) 

d 

Note that in the regime (|25| ) the subgap conductance decreases with the decrease of sheet conductance g. Comparing 
it with Eq.(||) we conclude that Gyi(ff) passes via maximum as g decreases. Using Eq. (|To|) together with Eq. (p5|) , we 
can find, at A — and sufficiently low frequencies, that phase fluctuations grow as 

S Q (t) = In 2 u d t (26) 

47T.9 

(note the similarity of Eq.(|2^) with the expression for the tunnelling action [^9| describing zero-biased anomaly of 
tunnelling conductance into a dirty 2D metal J|^]). Now let us come back to the discussion of b c (g) dependence. We 



expect that at low enough g critical distance b c will scale as \nb c /d cx yfg. This estimate |32| follows from Eq.(|26[) and 
is analogous to Finkelstein's result j^] g c ~ V^bcs = ^(^d/T®) for vanishing of T c in a homogenously disordered 
film. Let us now "turn on" repulsion constant A > 0, but assume that it is small, so the inequility opposite to Eq.([2^) 
is obeyed. As long as A is also less than the relevant scale <? -1 / 2 of the inverse logarithm in (|25|), the result valid 
for A = holds. However, there exists additional, intermediate case when g^ 1 ^ 2 < A < Gt/2irg. Here we expect 
the logarithm in ( ^5|) to be replaced by A" 1 , and therefore Ga ~ ffA. Thus we arrive at the following picture of the 
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IZ(g) = \a.b c (g)/d evolution: it grows with g decrease according to Eq.(pl|), untill g reaches g\ « Gt/2nX where TZ{g) 
passes via the maximum value of order Gt/4, and then decreases as TZ(g) ~ gX in the region A -2 < g < Gt/2irX. 
Finally, at lowest g < A~ 2 it decreases as 11(g) ~ y/g. 

In all cases considered, IZ(g) <C g, which means that weak localization corrections to the conductance of 2D 



metal 1 30 3l]] are indeed weak even at b > b c . Therefore, upon destruction of superconductivity with increase of 
inter-island separation, the normal metal ground-state will be produced. 

Now we turn to the discussion of the low-temperature phase boundary between superconductive and metal states 
(assuming that the condition (|2^) is fulfilled, and, also, the temperature T is so low that single-electron transport 
via NS barriers can be neglected). We calculate the dependence T c (b) of the critical temperature on the separation 
between islands, assuming it is much less than the zero-T critical value b c from Eq.(pl|). As a result, we obtain values 
of T c which are higher than the effective charging energy E* c defined in Eq. (|l4|) . On the other hand, we will always 
assume T c <C u>\. Then the MFA transition temperature is determined by the equation 



The strength of Josephson coupling depends on temperature (cf.Eq.(|lq)), and this dependence must be taken into 
account when calculating the sum over islands J*r, which replaced J from Eq.(|2(]) at T > 0. At lower enough 
temperatures this effect can be accounted for by simply choosing the upper limit of integration in Eq. ( p0| ) to be the 
thermal coherence length £t = y/D/T. This results in 

In 

Jt = j £ie ; 6<<Ct (28) 



Solving Eq.(g7|) together with Eq.(g8[) we find 

T MFA (b) - D (29) 
c W ~fo 2 Aln 2 ^ [ ' 

This formula can be applied as long as £t ^S> b, i.e., when T^ SF (b) <C D/b 2 , which is equivalent to the condition 
b ^> b\ with 

(30) 

At shorter b, transition temperature grows in such a way that £t becomes shorter than 6, and the MFA equation 
( p7j ) cannot be used anymore. In this case an appropriate approximation will be to neglect all couplings between 
non-nearest neighbours, i.e. to consider Berezinsky-Kosterlitz-Thouless (BKT) phase transition in 2D XY model with 
interaction between nearest neighbours |33]| . Moreover, in this region quantum phase fluctuations becomes weak and 
can be neglected in determination of T c . Therefore we estimate transition temperature as T c « J(b,T c ) 1 where the 
nearest-neighbours Josephson coupling J(b, T) in the relevant temperature range can be evaluated as 



T(hT\~ G * 1 exp(-6/g T ) , h f r „ n 
J(6 ' T) ~8^^ A 2 ln 2 «f ' 6>>CT (3) 



Using Eq.(^l|) and solving the equation T c = J(b, T c ) for T c with logarithmic accuracy, we finally obtain 




The expressions (^9|) and ( p2[ ) approximately match each other at b ps b\ defined by Eq.flBOp. The phase transition is 
of mean-field kind at b 3> b\ whereas it is of the BKT type |33| at b <C b\ . Therefore the temperature dependencies 
of macroscopic superconductive parameters such as superfluid density p s {T) and critical current density j c (T) are 
expected to differ significantly depending on whether b is small or large with respect to b\. 

4. Conclusions. We have shown that in a mixed superconductor-normal system with a small fraction of super- 
conductive regions quantum fluctuations effects still persist, and lead to the destruction of superconductivity even at 
T = 0, in spite of the good metal sheet conductance a ^> 4e 2 /h and the low SC-M interface resitance R t <§C /i/4e 2 . 
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These effects decay exponentially with the Andreev subgap conductance Ga = G 2 /2TrgX 3> 1; in particular, macro- 
scopic superconductivity persists down to the relative fraction of SC material n sc w {d/b c ) 2 ~ exp(— G^/2 — 1/A) . 
The temperature scale relevant for the effect we discussed is about T\ = D/b 2 .- At the distances between neighboring 
islands b <C b c the SC-M transition line T c (b,g) is determined by Eq.(p9f) and Eq.(|3^). These results were derived 
under the assumption that the Andreev subgap resistance is mainly determined by the S-N boundary. With the 
decrease of g, the subgap resistance first decreases but then starts to grow when the resistance of the film dominates 
the S-N boundary resistance. Quantitative investigation of quantum fluctuations and phase diagram in the corre- 
sponding region of parameters (the condition opposite to Eq.(^)) will be a subject of future work; we expect that 
at smallest g < A~ 2 the universal result — lnn sc cx yfg holds. It is important that the considered effects of quantum 
phase fluctuations are strong in the temperature range there usual interference and interaction-induced corrections 
to conductivity are still weak. In the normal metal state and at relatively high temperatures these fluctuations are 
expected to produce logarithmic in temperature corrections to conductivity, with the magnitude considerably larger 
than the weak-localization ones. 

There are several open problems closely related to the simple model we have studied here. It would be interesting 
to study the influence of quantum (T=0) fluctuations of the phases <frj(t) in the metallic phase ( b > b c ) on weak 
localization corrections to metal conductance and to check if they produce finite "dephasing time" t^(0). Strong 
dependence of the effective Coulomb energy Eq on the film conductance indicate an importance of mesoscopic fluctu- 
ations of conductance between different regions of nominally uniform film. Intersting problem is related to the effect 
of magnetic field of the scale B ~ $o/b 2 , which may lead to a formation of a superconductive glass state (cf. e.g. fjijl). 
Another important problem is the magnetic field - induced SC-normal transition in a uniformly disordered SC film at 
low T: it was argued by Spivak and Zhou Q , that mesoscopic fluctuations leads to a finite probability to find small 
SC cluster under the field which is much above "mean" H C 2- The problem of phase coherence between such clusters 
would resemble our model under the external magnetic field, but with important modifications due to randomness in 
the cluster positions. 
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